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Genistein and bisphenol A exposure cause estrogen 
receptor 1 to bind thousands of sites in a cell 
type-specific manner 
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Endogenous estrogens that are synthesized in the body impact gene regulation by activating estrogen receptors in diverse 
cell types. Exogenous compounds that have estrogenic properties can also be found circulating in the blood in both 
children and adults. The genome-wide impact of these environmental estrogens on gene regulation is unclear. To obtain 
an integrated view of gene regulation in response to environmental and endogenous estrogens on a genome-wide scale, we 
performed ChlP-seq to identify estrogen receptor 1 [ESR1; previously estrogen receptor ol) binding sites, and RNA-seq in 
endometrial cancer cells exposed to bisphenol A (BPA; found in plastics), genistein [GEN; found in soybean), or 17p- 
estradiol (E2; an endogenous estrogen). GEN and BPA treatment induces thousands of ESR1 binding sites and >50 gene 
expression changes, representing a subset of E2-induced gene regulation changes. Genes affected by E2 were highly 
enriched for ribosome-associated proteins; however, GEN and BPA failed to regulate most ribosome-associated proteins 
and instead enriched for transporters of carboxylic acids. Treatment-dependent changes in gene expression were associ- 
ated with treatment-dependent ESR1 binding sites, with the exception that many genes up-regulated by E2 harbored a BPA- 
induced ESR1 binding site but failed to show any expression change after BPA treatment. GEN and BPA exhibited a similar 
relationship to E2 in the breast cancer line T-47D, where cell type specificity played a much larger role than treatment 
specificity. Overall, both environmental estrogens clearly regulate gene expression through ESR1 on a genome-wide scale, 
although with lower potency resulting in less ESR1 binding sites and less gene expression changes compared to the en- 
dogenous estrogen, E2. 



[Supplemental material is available for this article.] 

Endogenous estrogens are physiologically important signaling 
molecules that can vary in concentration by an order of magnitude 
within a reproductive female. Despite such large fluctuations, es- 
trogens play pivotal and diverse roles in many tissues. In addition 
to endogenous estrogens, foreign estrogenic compounds can be 
found in the environment as well as at high doses in many people's 
blood streams (Grace et al. 2004; Jeng et al. 2010), providing the 
opportunity for these exogenous compounds to alter estrogen 
signaling in a variety of cell types. 

Epidemiological studies suggest that exposure to environmen- 
tal estrogens affects cancer rates (Ardies and Dees 1998) and impacts 
the age of puberty onset in young women (Aksglaede et al. 2009). 
Two compounds of particular interest due to their abundance are 
bisphenol A (BPA), used in the construction of polycarbonate plas- 
tics as well as epoxy resins found in water bottles, food containers, 
and dental fillings, and genistein (GEN), a prevalent isoflavone 
found in soybean. Both GEN and BPA can be observed at high levels 
(>40 |jlM) in circulating serum (Holder et al. 1999; Schonf elder 
et al. 2002), suggesting that these compounds have the oppor- 
tunity to alter estrogen signaling and gene regulation through 
estrogen receptors. 
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Observations in animal models have shown that both GEN 
and BPA can have a direct impact on development and cancer 
progression. Prenatal exposure to both GEN and BPA causes an 
acceleration in puberty onset (Nikaido et al. 2004) and desensitizes 
tissues in adult animals to the endogenous estrogen, 17(3-estradiol 
(E2) (Woodruff and Walker 2008; Molzberger et al. 2011). BPA 
exposure increases the risk of breast cancer in animal models 
(Weber Lozada and Keri 2011), while dietary GEN is able to block 
tamoxifen-induced growth inhibition of breast tumors (Du 
et al. 2012). 

While there are clear physiological effects of xenoestrogen 
exposure, an integrated molecular picture detailing how environ- 
mental estrogens act in combination with estrogen receptor 1 
(ESR1; previously estrogen receptor a) to impact gene regulation 
remains to be uncovered. GEN, along with some other environ- 
mental estrogens, exhibits binding affinities for purified ESR1 that 
are similar to E2 (Kuiper et al. 1998) and has been shown to induce 
ESR1 binding at particular transcriptional promoters (Chang et al. 
2008). BPA exposure has been shown to induce some of the gene 
expression changes that are also observed following treatment 
with E2 (Hess-Wilson et al. 2006, 2007), while eliciting some dis- 
tinct gene expression changes (Singleton et al. 2006). To appreciate 
the full impact of these exogenous chemicals on gene regulation, 
an integrated approach that includes the investigation of genome- 
wide changes in ESR1 binding and gene expression caused by en- 
vironmental estrogen exposure in multiple cell types is needed. 

To determine the gene regulatory consequences of treatment 
with exogenous and endogenous estrogens, we exposed endometrial 
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cancer cells to GEN ; BPA, or E2 and mea- 
sured ESR1 binding and gene expression 
changes. We found that GEN and BPA 
cause a distinct set of transcriptional re- 
sponses that are observed after E2 treat- 
ment. Similar patterns were also observed 
in breast cancer cells, which exhibited 
very different gene regulatory responses to 
estrogens but maintained the same re- 
lationship between E2, GEN, and BPA. 
Importantly BPA induced a small number 
of gene expression changes, yet influenced 
E2-dependent gene expression, suggesting 
a mechanism whereby endogenous estro- 
gen signaling is altered in the presence of 
environmental estrogens. 

Results 
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E2-induced ESR1 Binding Site Enrichment 



BPA and GEN promote ESRl binding 
to more than 1000 genomic loci 

The genomic pathway of estrogen signaling consists of estrogen re- 
ceptors binding to genomic loci and regulating transcription. To 
determine the impact that environmental estrogens have on the 
genomic pathway we used ChlP-seq to identify sites of occupancy 
for ESRl following treatment with GEN and BPA for 1 h. We also 
measured binding of ESRl after 1 h of E2 treatment to compare the 
effects of environmental estrogens to those of an endogenous es- 
trogen. Using 1-h DMSO treatment followed by ChlP-seq targeting 
ESRl as a vehicle control, we used MACS (Zhang et al. 2008) to 
discover ESRl binding sites induced by each compound in two 
independent biological replicates of the endometrial cancer cell 
line, ECC-1, which proliferates in response to E2 exposure (Castro- 
Rivera et al. 1999). For E2, GEN, and BPA, only sites that were 
replicated in each experiment were analyzed further. 

There were reproducible differences in the number of ESRl 
binding sites detected after treatment with each chemical. There 
were 8622 ESRl -bound sites induced by E2, half (4910) as many 
sites induced by GEN, and one-fifth (1471) as many sites induced 
by BPA (Table 1). This is the first evidence that BPA and GEN ex- 
posure causes ESRl to bind thousands of sites across the genome. 
To determine, relative to E2, whether GEN and BPA induction 
caused ESRl to bind to distinct sites, we analyzed the overlap of 
ESRl binding sites associated with each compound. Almost all of 
the sites bound after GEN and BPA treatment were also bound after 
E2 treatment, with 93.6% of GEN-induced sites and 95.2% of BPA- 
induced sites overlapping with E2-induced sites. 

Figure 1, A and B, shows two prototypical sites in which ESRl 
binding is induced by all three ligands. The overall ESRl binding 



Table 1. Number of ESRl binding sites and overlap 



Figure 1. ESRl binding is induced by GEN and BPA. Examples of ESRl binding induced by E2, GEN, 
BPA, and DMSO near the promoters of GREB1 (A) and JUND (B) in ECC-1 are shown. The height of the 
orange graph represents number of reads per 1 million total reads. (C) Binding site signal strength is 
plotted for each binding site with E2-induced signal strength (replicate 1 ) on the x-axis vs. E2 (black; 
replicate 2), GEN (green; replicate 1), and BPA (red; replicate 1) signal strength on the /-axis for 
ECC-1. 
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signal after E2 induction is greater than the signal after GEN in- 
duction, which is greater than the signal after BPA induction. To 
examine the quantitative differences in ESRl binding induced by 
different compounds in more detail, we compared ESRl binding 
signal strength, calculated as the number of sequence reads in 
a binding site per million aligned reads (RPM), for every binding 
site in each ChIP experiment. For example, the signal strengths for 
E2, GEN, BPA, and DMSO of the binding site shown in Figure 1A 
are 21.5, 13.2, 4.1, and 0.15 RPM. Figure 1C displays the signal 
strengths for all E2-induced ESRl binding sites. We found that the 
strength of ESRl binding after E2 treatment was very similar be- 
tween biological replicates (r > 0.85). We also found that the 
strongest signal binding sites after E2 induction had the strongest 
signal after GEN and BPA induction (r > 0.8), with GEN-induced 
sites averaging 1.66-fold lower signal and BPA-induced sites aver- 
aging 5.1 -fold lower signal compared to E2-induced sites. 

To identify differences in ESRl binding sites induced by each 
compound, we split ESRl binding sites into three groups: sites 
induced by E2 only, sites induced by E2 and GEN (E2+GEN), and 
sites induced by all three compounds (E2+GEN+BPA). The distri- 
bution of ESRl binding sites relative to gene features is nearly 
identical for each group (Fig. 2A). The majority of ESRl binding 
sites are found in intergenic regions at least 2 kb away from known 
genes. When compared to the genomic background, each group of 
binding sites was enriched for regions near transcription start sites: 
promoters, first exons, and first introns (P < 1 X 10~ 6 ; x 2 test). A 
common sequence motif was found in all three binding site groups, 
which matches the canonical estrogen response element (ERE) 
(Fig. 2B). We analyzed each binding site for the best match to the motif 
and estimated a relative affinity of ESRl for the site (Hertz and Stormo 
1999). Figure 2C shows that each set of binding sites is enriched for 
strong matches to the ESRl motif compared to randomly chosen se- 
quences. ESRl sites bound after all three treatments as well as sites 
bound after GEN and E2 treatment have significantly higher scoring 
matches than sites induced by E2 only (P < 10~ 4 ; Wilcoxon test). 
Overall, ESRl binding sites induced by each environmental es- 
trogen represent a subset of the sites induced by the endogenous 
estrogen, E2, with GEN- and BPA-induced sites exhibiting an 
enrichment for strong matches to an ERE and higher signal after 
E2 treatment. 
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Figure 2. Properties of E2-, GEN-, and BPA-induced ESR1 binding sites. (A) The proportion of genomic locations of ESR1 binding sites with respect to 
genes is displayed. The genomic background is also shown. (S)The most prevalent sequence motif found in ESR1 binding sites found after E2-, GEN-, or 
BPA-treatment matches the canonical estrogen response element (ERE). (C) Distributions of ERE motif strength, measured as log relative affinity 



compared to highest possible affinity site, are shown as boxplots. Boxes cover 25 th 
interval. 



to 75 percentiles, and the whiskers define the 95% confidence 



E2, GEN, and BPA cause overlapping gene expression responses 

Given that both BPA and GEN cause ESR1 to bind many sites across 
the genome, we sought to determine changes in gene expression 
that each environmental estrogen induced. We treated ECC-1 with 
E2, GEN, BPA, or DMSO (as a control) for 8 h, a time point at which 
an even representation of ESR1 binding sites at up- and down- 
regulated genes has been observed (Carroll et al. 2006). We then 
performed RNA-seq using the Tn-RNA-seq protocol that we recently 
developed (Gertz et al. 2012), obtaining more than 40 million 
paired-end 50-bp reads for each sample. Biological replicates were 
performed for each treatment in order to determine significantly 
affected transcripts (see Methods). 

In total, more than 900 genes were significantly affected by E2 
treatment, with a false discovery rate of 5% (see Supplemental 
Tables SI and S2 for gene expression data). In response to E2 ex- 
posure, 547 genes were up-regulated and 440 genes were down 
regulated. This includes the down-regulation of one long inter- 
genic noncoding (LINC) RNA, LINC00152, as well as up-regulation 
of one microRNA precursor, miR-663a, thought to play a role in 
nasopharyngeal carcinoma through the repression of CDNK1A (Yi 
et al. 2012). Genes whose expression is affected by E2 treatment are 
highly enriched in ribosomal-associated proteins that play a role in 
translation (Supplemental Table S3), which is consistent with the 
proliferative effect of E2 on ECC-1. 

Congruent with the observation that GEN and BPA induce 
fewer ESR1 binding sites than E2, both compounds affect a smaller 
number of genes in comparison to E2. Treatment with GEN caused 
up-regulation of 207 genes and down-regulation of 49 genes. 
While 54 genes were up-regulated after BPA exposure, only two 
genes were down-regulated. Unlike transcripts that were affected 
by E2, genes induced by GEN and BPA are not enriched for 
translation and ribosome function. In both cases, genes in- 
volved in transport of carboxylic acids were overrepresented (P < 
2.51 X 10~ 5 ). This includes the transport of amino acids, sialic acid, 
prostaglandins, and other mono-, di-, and tricarboxylic acids. 



The genes whose expression was influenced by GEN or BPA 
exposure represent a subset of the genes impacted by E2 treatment. 
Of the genes that were up-regulated by GEN and BPA, 94.2% and 
94.4% were also up-regulated by E2, respectively (Fig. 3 A). GEN and 
BPA down-regulated genes overlap with E2 down-regulated genes at 
a rate of 80% and 100%, respectively. There were no genes that are 
affected by E2 and regulated in the opposite direction by GEN or 
BPA. The high overlap between genes affected by E2, GEN, and BPA 
indicates there is some functional redundancy between these com- 
pounds, with the environmental estrogens acting on a subset of 
genes influenced by E2 that is depleted for ribosomal function and 
enriched for carboxylic acid transport. A small set of genes is sig- 
nificantly impacted by GEN or BPA and not E2; however, each of 
these genes responds to E2 in the same manner as GEN or BPA at 
levels that are just below significance. 

BPA and GEN also induce smaller magnitude gene expression 
changes in comparison to E2. On average, gene expression changes 
brought on by GEN are 10% smaller than changes caused by E2 
treatment (P < 2.2 X 10~ 16 ; Wilcoxon test), and expression changes 
induced by BPA are 25% smaller than E2 changes (P < 2.2 X 10" 16 ; 
Wilcoxon test) (Fig. 3B). Figure 3C shows some prototypical ex- 
amples of genes significantly affected by E2. TGFA, a gene known 
to be regulated by E2 (Fontana et al. 1992) and overexpressed in 
endometrial cancer (Niikura et al. 1996), is up-regulated by E2, 
GEN, and BPA. RRP1B, up-regulated only by E2, is thought to play 
a role in ribosomal RNA processing and may be involved in met- 
astatic potential (Crawford et al. 2007). SLC22AS, a carnitine 
transporter associated with autoimmune diseases (Peltekova et al. 
2004; Santiago et al. 2006), is significantly up-regulated by E2 only 
but shows similar changes after GEN and BPA treatment. IGFBP3, 
a modulator of insulin-like growth factor (Firth et al. 1998) and 
a direct target of androgen receptor (Peng et al. 2006), is significantly 
down-regulated by E2 and GEN; however, BPA down-regulates 
IGFBP3 to a level that does not meet genome-wide significance. 
The lower-signal ESR1 binding sites and lower-magnitude gene 
expression changes induced by BPA and GEN, compared to E2, 
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Figure 3. Gene expression changes in response to GEN and BPA. (A) The overlap between genes 
up-regulated (left panel) and genes down-regulated (right panel) by E2, GEN, and BPA. (B) A 
scatter plot shows the relative magnitude of gene expression changes with log 2 of the fold change 
induced by E2 on the x-axis, and log 2 of the fold change induced by GEN (green) and BPA (red) on 
the /-axis. The black line represents equal magnitude. (C) Examples of gene expression mea- 
surements in DMSO, E2, GEN, and BPA are shown for four representative genes. Error bars rep- 
resent s.e.m. 



suggest that these environmental estrogens 
are less potent inducers of gene regulation 
through ESR1. 

Many genes bound after BPA or GEN 
treatment are regulated by E2 but not 
by BPA or GEN 

Based on the overlap of both ESR1 bind- 
ing sites and gene expression changes 
induced by E2, GEN, and BPA, we sought to 
quantify the relationship between ESR1 
binding and gene expression changes after 
exposure to environmental estrogens. We 
partitioned both genes and ESR1 -bound 
sites into three groups: E2 only, E2+GEN, 
and E2+GEN+BPA. A breakdown of the 
number of genes and binding sites can be 
found in Supplemental Table S4. We ana- 
lyzed the genes in each group and whether 
they had an associated ESR1 binding site 
(within 50 kb of the transcription start 
site). Figure 4 shows the number of genes 
in each group that can be explained by 
nearby ESR1 binding sites. As expected, 
the highest association was found when 
the group of genes and the group of bind- 
ing sites matched. For example, 63% 
of E2+GEN+BPA genes had an associated 
E2+GEN+BPA binding site, compared to 
7.5% of E2+GEN+BPA genes being associ- 
ated with E2 only binding sites or 11% be- 
ing associated with E2+GEN binding sites. 

There are a large number of genes 
associated with a BPA-induced ESR1 bind- 
ing site (E2+GEN+BPA) which are not reg- 
ulated by BPA treatment (E2 only or 
E2+GEN). In total, 43 genes up-regulated 
by E2, and not by BPA, had an associated 
E2+GEN+BPA ESR1 binding site. MAP4K3, 
a nutrient sensitive regulator of mTOR 
signaling (Findlay et al. 2007), is up-regulated 
by E2 and not by BPA; however, BPA 
treatment as well as E2 treatment causes 
ESR1 to bind a site 20 kb downstream from 
the transcription start site (Fig. 4C). GEN 
treatment can also induce ESR1 binding 
without causing a significant gene expres- 
sion response, as 53 genes are associated 
with an E2+GEN or E2+GEN+BPA ESR1 
binding site without GEN affecting their 
gene expression. SKIL, a negative regu- 
lator of TGFp signaling (Lamouille and 
Derynck 2009), harbors three GEN-induced 
ESR1 sites within 20 kb of the transcrip- 
tion start site but does not respond to 
GEN treatment (Fig. 4C). The inability of 
BPA and GEN to produce gene expression 
responses, while causing ESR1 to bind 
nearby sites, suggests that BPA- and GEN- 
induced ESR1 binding sites are less potent 
mediators of gene regulation. 
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Figure 4. Association between ESR1 binding and gene expression changes in response to E2, GEN, 
and BPA. The number of genes in each group that are within 50 kb of ESR1 binding sites (also broken 
out by group) for up-regulated (A) and down-regulated (B) genes. (C) ChlP-seq signal and gene ex- 
pression levels after E2-, GEN-, BPA-, and DMSO-treatment are shown for two representative genes. The 
height of the orange graph represents the normalized number of reads for each ChlP-seq experiment. 
Error bars represent s.e.m. 



Gene expression changes are modulated by continuous 
exposure to BPA or E2 

To explore possible interactions between exogenous and endoge- 
nous estrogens, we focused on how E2 and BPA influence gene 
regulation in combination. We grew ECC-1 with continual expo- 
sure to either 100 nM BPA, 5 pM E2, 10 nM E2, or DMSO for 5 d. We 
then performed short-term treatment of cells with BPA or E2 and 
measured relative changes in gene expression for a representative 
panel of genes (Fig. 5). Continual treatment with 10 nM E2 is meant 
to mimic adult levels of estrogen, and the continual treatment with 
5 pM is meant to mimic prepubertal levels of estrogen (Klein et al. 
1999). 

We found that continuous exposure to E2 influences the gene 
expression response induced by BPA. Figure 5A shows that the 



expression of TGFA, ZNF608, DOK7, and 
SLC22A5 are all significantly affected by 
short-term exposure to BPA when cells are 
grown continuously with only DMSO. 
The expression changes observed in this 
assay are very similar to the results obtained 
with RNA-seq (r = 0.85). When ECC-1 
cells are continuously exposed to pre- 
pubertal levels of E2, BPA-induced expres- 
sion changes remain significant. However, 
when cells are pretreated with adult levels 
of E2, short-term BPA exposure does not 
change gene expression. These results 
indicate that chronic adult E2 exposure 
blunts the BPA-induced gene expression 
changes. By contrast, preadolescent levels 
of E2 do not diminish gene expression 
changes mediated by acute BPA treatment. 

Interestingly, in the context of chronic 
exposure to BPA, we found attenuation in 
the gene expression changes induced by 
E2. Adult levels of E2 were able to induce 
significant gene expression changes for 
TGFA, ZNF608, DOK7, and SLC22A5 in 
cells with long-term exposure to pread- 
olescent levels of E2, whereas in cells pretreated with BPA, E2 did 
not elicit a significant change in the expression of TGFA, ZNF608, 
and SLC22A5 (Fig. 5B). The fold change of DOK7 expression in 
response to adult levels of E2 was reduced by half with continuous 
exposure to BPA. Long-term exposure to BPA enhanced basal E2- 
indpendent gene expression of TGFA, ZNF608, SLC22AS, and to 
a lesser extent, DOK7, to levels similar to those observed upon E2 
treatment. This suggests that long-term BPA exposure influences 
E2-dependent gene expression. 

To examine whether the impact of long-term exposure to BPA 
on gene expression is ESR1 -dependent, we treated ECC-1 cells with 
BPA and the ESR1 antagonist ICI 182,780 (Fulvestrant) and mea- 
sured the expression of TGFA, ZNF608, SLC22A5, and DOK7 
(Supplemental Fig. SI). We found that for each gene, the expres- 
sion changes driven by long-term BPA exposure are lost when cells 
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Figure 5. Long-term exposure to BPA and E2 alters gene expression response. (A) Relative expression measurements are shown for four genes in ECC-1 
after short-term BPA exposure following 5 d of treatment with DMSO, 5 pM E2, or 10 nM E2. (B) Relative expression measurements are shown for four 
genes in ECC-1 after short-term treatment with 10 nM E2 following 5 d of treatment with DMSO, 5 pM E2, or 5 pM E2 plus 100 nM BPA. Error bars 
represent s.e.m. 
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are treated with ICI. These results indicate that long-term BPA 
treatment acts through ESR1 to influence gene expression. 

BPA and GEN exhibit similar gene regulation patterns in T-47D 

To determine the impact of environmental estrogens on ESR1 
signaling in another cell type, we performed ChlP-seq and RNA- 
seq after DMSO, E2, GEN, or BPA treatment in the breast cancer cell 
line T-47D. We found a striking amount of cell type specificity in 
ESR1 binding and gene expression changes induced by E2, GEN, 
and BPA. While both cell lines had a similar number of ESR1 
binding sites, 7453 in T-47D and 8622 in ECC-1, only 1446 of the 
ESRl-bound sites overlap (19% of T-47D bound sites) (Fig. 6A). 
These results indicate that the set of ESRl-bound sites is distinct 
between endometrial and breast cancer cells, consistent with pre- 
vious reports indicating high levels of cell type specificity in ESR1 
binding (Krum et al. 2008; Hurtado et al. 2011). 

The disparity in ESR1 binding could be due to differences 
between breast and endometrial cells specifically or to differences 
between cells regardless of type. To address this question, we com- 
pared, by ChlP-seq, the sites bound by ESR1 identified in ECC-1 and 
T-47D to ESR1 binding sites found in MCF-7, a breast cancer cell 
line (Welboren et al. 2009). Of the 10,196 ESR1 binding sites 
identified in MCF-7 cells, T-47D shares 3037 (41% of T-47D bind- 
ing sites) and ECC-1 shares 1798 (21% of ECC-1 binding sites). The 
two breast cancer cell lines, MCF-7 and T-47D, are twice as similar 



to each other than they are to the endometrial cancer line, in- 
dicating that at least these two breast cancer cell lines share some 
common genomic traits that result in ESR1 binding to similar lo- 
cations (Supplemental Fig. S2). While the breast cancer cell lines 
are more similar to each other, there remain 4022 (53.9%) ESR1 
binding sites specific to T-47D, suggesting that cell-autonomous 
regulation of ESR1 binding also occurs. 

Very similar relationships between ESR1 binding after treat- 
ment with E2, GEN, and BPA were observed in ECC-1 and T-47D 
(Table 1). E2 induced twice as many ESR1 binding sites as GEN and 
five times as many sites as BPA in T-47D. Binding sites induced by 
GEN and BPA overlapped highly (>91%) with E2 induced sites, 
indicating that BPA and GEN treatment causes binding of ESR1 to 
a subset of the sites bound following E2 treatment. In fact, GEN- 
and BPA-induced ESR1 binding matched much more closely to E2- 
induced sites in the same cell line than GEN- or BPA-induced sites 
in the other cell line (—90% vs. —20%). Consistent with data from 
ECC-1, ESR1 binding site signal strengths in T-47D were strongest 
after E2 treatment, lower after GEN treatment, and lowest after BPA 
treatment (Fig. 6B). In both ECC-1 and T-47D, GEN and BPA in- 
duced a smaller number of ESR1 binding sites, compared to E2 but 
maintained a high level of cell type specificity. 

FOXA1 and GATA3 are transcription factors that have been 
shown to be important for ESR1 binding and enhancer activity in 
breast cancer cells (Kong et al. 2011). To determine if ESR1 binding 
sites induced by E2, GEN, and BPA were more or less likely to co-occur 
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Figure 6. GEN and BPA induce ESR1 binding and expression changes in T-47D. (A) The overlap between E2-induced ESR1 binding sites in ECC-1 and 
T-47D is displayed. (B) Binding site signal strength is plotted for each binding site with E2-induced signal strength (replicate 1 ) on the x-axis vs. E2 (black; 
replicate 2), GEN (green; replicate 1), and BPA (red; replicate 1) signal strength on the /-axis for T-47D. (C)The overlap between genes regulated by E2, 
GEN, or BPA is shown as a Venn diagram. (D) The number of genes in each group that are within 50 kb of ESR1 binding sites (also broken out by group) is 
shown for up-regulated genes. 
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with these factors, we mapped FOXA1 and GATA3 binding by 
ChlP-seq in untreated (DMSO control) T-47D cells. We found a sig- 
nificantly higher percentage of E2+GEN+BPA (39%) or E2+GEN 
(32.7%) ESR1 sites co-occuning with FOXA1 and GATA3 compared 
to E2-only ESR1 sites (26%) (P < 2.2 X 1(T 16 and P < 2.7 X 1(T 10 , 
respectively; Fisher's exact test). The higher overlap between ESR1 
sites, found after BPA and GEN treatment, and FOXA1 and GATA3 
sites suggest that binding of ESR1 induced by GEN and BPA may 
have higher dependency on cof actors. 

In a pattern similar to ECC-1, GEN and BPA caused a subset of 
the gene expression changes induced by E2 in T-47D. E2 treatment 
led to the up-regulation of 443 genes and the down-regulation of 
168 genes, while GEN up-regulated 106 genes and down-regulated 
14 genes (Fig. 6C). BPA only affected one gene, up-regulating 
CXCL12. CXCL12, a chemokine that can increase metastatic po- 
tential (Muller et al. 2001), was up-regulated by E2 and GEN as 
well. We observed high overlap between GEN- and E2-affected 
genes; every gene that was up-regulated by GEN in T-47D was also 
up-regulated by E2, and 10 out of 14 genes down-regulated after 
GEN treatment were also down-regulated by E2. Genes affected by 
both GEN and E2 were enriched for tissue development (P < 7.53 X 
10~ 4 ), and genes affected by E2 were enriched in ribosomal pro- 
cesses and translation (P < 2.50 X 10~ 24 ). This finding is consistent 
with ECC-1 and the increase in proliferation that has been ob- 
served after E2 treatment in T-47D (Chalbos et al. 1982). 

While E2 regulates ribosomal genes in both ECC-1 and T-47D, 
the overlap between genes regulated by E2 in each cell line is small. 
Consistent with ESR1 binding, gene expression changes induced 
by each compound were highly cell type-specific. Only 71 up- 
regulated genes (16% of T-47D up-regulated genes) overlapped 
between ECC-1 and T-47D; nine genes were down-regulated in 
both lines (5% of T-47D down-regulated genes). There were 46 
genes that changed direction between cell lines, with 39 genes up- 
regulated in T-47D while being down-regulated in ECC-1, and nine 
up-regulated in ECC-1 while being down-regulated in T-47D. Of 
the genes regulated in opposite directions, 10 are ribosomal pro- 
teins (P < 5.8 X 10~ 17 ); however, most ribosomal proteins are up- 
regulated in one or both lines. 

We analyzed the overlap between gene expression changes 
and ESR1 binding in T-47D using the same groups described above: 
E2 only, E2+GEN, and E2+GEN+BPA. We found very similar pat- 
terns in T-47D as was observed in ECC-1 (Supplemental Table S4), 
with each group of genes being enriched for the same group of 
ESR1 binding sites (Fig. 6D). There were 72 genes that were asso- 
ciated with a BPA-induced ESR1 binding site which showed no 
expression change after treatment with BPA and a significant ex- 
pression change after E2 exposure. GEN caused ESR1 binding at 46 
E2-affected genes that did not change after GEN treatment. These 
results are consistent with the hypothesis that BPA- and GEN- 
induced ESR1 binding events are less potent mediators of gene 
regulation. 

Discussion 

By studying exposure of ECC-1 to GEN and BPA, we were able to 
demonstrate that these compounds induce thousands of ESR1 
binding sites. To our knowledge, this is the first genome-wide study 
of ESR1 binding in response to BPA or GEN treatment. We also 
found that GEN and BPA change the expression of a subset of genes 
affected by E2, representing 26% and 6%, respectively. Genes af- 
fected by GEN and BPA were more often up-regulated, as <20% 
were down-regulated. The tendency of GEN and BPA to cause 



up-regulation as opposed to down-regulation could be the result of 
different mechanisms. It is possible that ESR1 is more likely to bind 
coactivators than corepressors when induced by GEN or BPA. GEN 
and BPA may also be working on a longer time scale, since the 
majority of early gene expression changes induced by E2 involve 
up-regulation (Frasor et al. 2003), the gene regulation response 
caused by GEN and BPA after 8 h may be similar to early changes 
brought on by E2. 

There are many genes that are regulated by E2 only but have 
nearby ESR1 binding induced by E2, GEN, and BPA. This suggests 
that environmental estrogen-induced ESR1 is a less potent medi- 
ator of gene regulation, possibly a consequence of proportionally 
less ESR1 binding caused by either lower amounts of activated 
ESR1 or lower affinity ESR1 binding. The failure of GEN and BPA to 
cause significant gene expression responses could also be reflective 
of a failure to recruit coactivators. This hypothesis is consistent 
with in vitro data showing that BPA- and GEN-liganded ESR1 have 
100,000- and 10,000-fold lower relative recruitment of NCOA1 
and NCOA2 compared to E2, respectively (Routledge et al. 2000). 
The lack of recruitment of coactivators to BPA- and GEN-liganded 
ESR1 may explain the disconnect between ESR1 binding observa- 
tions and gene expression changes. Overall, there are many more 
ESR1 binding sites compared to gene expression changes after each 
treatment, in agreement with previous studies (Cheng et al. 2009; 
Reddy et al. 2009; Simpson et al. 2012). It is unclear if the extra 
binding sites represent spurious binding that is unrelated to gene 
expression or if most genes require multiple long distance binding 
sites. 

Genes impacted by E2 were highly enriched for ribosome- 
associated proteins, consistent with the observation that estrogen 
exposure increases uterine ribosome concentration in vivo (Moore 
and Hamilton 1964). However, GEN and BPA were unable to reg- 
ulate ribosome-associated genes. This raises the question of the 
influence that each compound has on proliferation, as making 
ribosomal proteins is an important step in increasing proliferation. 
In agreement with this observation, BPA has been reported to have 
no impact on ECC-1 proliferation at the dose used in our study 
(Bergeron et al. 1999). Previous studies of GEN's impact on pro- 
liferation suggest that intermediate doses, such as the 100 nM 
concentration that we used, increase proliferation (Hirsch et al. 
2007). GEN's failure to up-regulate ribosomal proteins may be due 
to the time point that we assayed or a difference in the mechanism 
with which GEN induces proliferation compared to E2. For example, 
GEN induces the expression of cyclin D2 and API component FOS, 
which are not affected by BPA and could be responsible for in- 
creased cell proliferation. 

BPA and E2 seem to desensitize ECC-1 to each other. When 
BPA was introduced to cells that were already exposed to high 
levels of E2, there was no effect on gene expression. However, 
when cells were pretreated with low amounts of E2, similar to 
prepubertal levels, BPA was able to influence gene expression. Our 
results suggest that BPA is affecting E2-dependent gene expression 
in the context of low E2 levels. This is consistent with findings that 
infants and children, who exhibit low concentrations of E2, are 
more sensitive to BPA exposure than are adults (Crain et al. 2008; 
Woodruff and Walker 2008). Long-term treatment with BPA is also 
able to block gene expression changes brought on by the switch 
from low to high E2 levels, suggesting that the presence of BPA may 
augment the molecular response to natural changes in estrogen 
concentrations. 

We also studied the influence of BPA and GEN on gene reg- 
ulation in a breast cancer cell line, T-47D. In comparing the cell 
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types, we found that E2 caused very different molecular responses 
in ECC-1 and T-47D, with <20% overlap of ESR1 binding sites and 
the genes impacted by E2. This is a surprising finding, as E2 induces 
proliferation in both cell lines (Chalbos et al. 1982; Castro-Rivera 
et al. 1999), and indicates that diverse gene regulatory responses 
can result in related phenotypic outcomes. While ECC-1 and 
T-47D differ substantially in their response to E2, the relationship 
between E2, GEN, and BPA within a line was consistent. Different 
cell types exhibited much lower overlap than different estrogens 
(—20% vs. —90%). The biggest difference in each compound's 
behavior between cell types was that BPA was only able to induce 
one significant gene expression change (CXCL12) in T-47D. This 
observation is surprising since BPA treatment causes 1877 ESR1 
binding sites in T-47D and could be due to a T-47D-specific factor 
that modulates gene expression changes from a BPA-liganded 
ESR1. The mechanism underlying this lack of gene expression re- 
sponse in a cell type-specific manner is unknown and requires 
further investigation. 

Our results show that physiological concentrations of GEN 
and BPA induce ESR1 binding but affect gene regulation to a dif- 
ferent extent. These differences are consistent between cell types 
and congruent with the opposite affects that GEN and BPA seem to 
have on human health. For example, prenatal exposure to BPA is 
associated with morphological abnormalities, reduced fertility, 
and higher rates of breast and prostate cancer (Maffini et al. 2006). 
In contrast to BPA, consumption of GEN is correlated with de- 
creased cancer risk (Adlercreutz 1995). However, early exposure to 
GEN is thought to delay the onset of puberty ( Jacobson-Dickman 
and Lee 2009) and lead to reduced fertility and underdevelopment 
of the female reproductive tract (Jefferson et al. 2007). Gene reg- 
ulatory differences between GEN and BPA may be sufficient to 
impart distinct physiological characteristics to the cell, or, alter- 
natively, work in combination with other factors to promote 
physiological differences upon exposure to these environmental 
estrogens. For example, although neither ECC-1 nor T-47D cells 
express estrogen receptor 2 (formerly estrogen receptor (3) under 
our culture conditions, estrogen receptor 2 has been shown to be 
growth-inhibitory, promote differentiation (Korach et al. 1996; 
Forster et al. 2002; Strom et al. 2004), and has higher affinity for 
GEN compared to ESR1 (Kuiper et al. 1998). This may explain some 
of the opposing roles of GEN and BPA in human health. GEN and 
BPA may also exert their effects through other signaling pathways. 
BPA has been shown to activate MAPK1 (also known as ERK) kinase 
(Lee et al. 2008), AKT1 (also known as AKT) kinase, and PIK3CA (also 
known as PI 3 -kinase) (Masuno et al. 2005) in different cell types, 
while GEN inhibits tyrosine-specihc kinases such as EGFR (Akiyama 
et al. 1987). Understanding how these other mechanisms interact 
with estrogen signaling pathways will be important in appreciating 
the impact that environmental estrogens have on human health. 

Methods 

Cell culture and inductions 

We grew both ECC-1 (ATCC CRL-2923) and T-47D (ATCC HTB- 
133) in RPMI-1640 (Invitrogen) supplemented with 1% penicillin- 
streptomycin (Invitrogen) and 10% fetal bovine serum (HyClone). 
Five days before induction, cells were switched to phenol red-free 
RPMI-1640 supplemented with 1% penicillin-streptomycin and 
10% Charcoal/Dextran-treated fetal bovine serum (HyClone). 
Fresh media was added 1 d before inductions. For ChlP-seq ex- 
periments, 6 (jlL of 5000 X concentrated solutions of E2, GEN, and 
BPA dissolved in DMSO were added to 30 mL media in 150-mm 



dishes. For the vehicle control, 6 \xL DMSO was added to 30 mL of 
media. For RNA-seq experiments, 3 |xL concentrated solution or 
DMSO was added to 15 mL of media in 100-mm dishes. For qPCR 
experiments, 10 jjlL 50x concentrated solution or DMSO was 
added to 500 |xL of media in a single well of a 24-well plate. ICI 
182,780 was used at a final concentration of 100 nM. 

For chromatin IP, chromatin was harvested by adding form- 
aldehyde to the media to a final concentration of 1% for 10 min. 
After 10 min, glycine was added to a final concentration of 0. 125 M 
for 5 min. Cells were washed with cold PBS, and then Farnham 
lysis buffer (5 mM PIPES at pH 8.0, 85 mM KC1, 0.5% NP-40) was 
added. Cell lysate was then collected with a policeman, and nuclei 
were collected by centrifugation for 5 min at 1000g. Pelleted nuclei 
were stored at -80°C. 

For RNA samples, a 100-mm dish of cells was lysed with 600 
jjlL Lysis/Binding buffer from an mRNA direct kit (Invitrogen). 
Lysate was then collected with a policeman and passed through 
a 20-gauge needle 10 times. Lysate was stored at -80°C. 

Chromatin immunoprecipitation and analysis 

Chromatin immunoprecipitation and ChlP-seq were performed as 
previously described (Reddy et al. 2009). The antibodies ERa (HC- 
20), GATA-3 (HG3-31), and HNF3a/0 (C-20) (FOXA1) (Santa Cruz 
Biotechnology, Inc.) were verified by immunoprecipitation fol- 
lowed by Western blot (data not shown) and were used in ChlP-seq 
experiments. At least 12 million 36-bp reads were obtained for 
every replicate, and each ChlP-seq experiment was performed in 
duplicate. 

To determine binding sites, sequences were aligned to the 
February 2009 (GRCh37/hgl9) build of the human genome using 
ELAND (Illumina). The alignments obtained from ESR1 ChlP-seq 
libraries of cells induced by E2, GEN, or BPA were compared to 
sequence alignments obtained from the DMSO induction ChlP- 
seq libraries in the same cell line using the program MACS (Zhang 
et al. 2008). FOXA1 and GATA3 ChlP-seq libraries were compared 
to reverse crosslink input control libraries from T-47D cells. At least 
78% of binding sites were shared between ChlP-seq replicates, and 
only binding sites that were identified in both replicates were used 
for further analysis. All ChlP-seq data can be accessed through the 
USCS Genome Browser's site "Transcription Factor Binding Sites 
by ChlP-seq from ENCODE/HAIB" (http://genome.ucsc.edu/cgi- 
bin/hgTrackUi?g=wgEncodeHaibTfbs) and through the Gene Ex- 
pression Omnibus (GEO) accession number GSE32465. 

RNA-seq and analysis 

Tn-RNA-seq was performed as previously described (Gertz et al. 
2012). In brief, mRNA was purified from cell lysate with an mRNA 
Direct kit (Invitrogen) according to the manufacturer's instructions. 
Purified mRNA was treated with 2 U DNase I (New England BioLabs) 
in RSB buffer (1 mM Tris pH 7.4, 1 mM NaCl, 300 nM MgCl 2 ) for 
10 min at 37°C. Two |xL of 500 mM EDTA pH 8.0 were added, 
and the reaction was purified using a spin-50 Sephadex column 
(USA Scientific). Double-stranded cDNA was constructed from the 
DNase-treated mRNA. Fifty ng cDNA then underwent Nextera 
tagmentation and PCR. Libraries were constructed with DNA 
barcodes and pooled as four per lane on an Illumina HiSeq 2000. 

We aligned paired-end 50-bp reads to spliced reference tran- 
scripts using bowtie (parameters: -n 2 -a -m 10 -X 3000). For each 
replicate, we obtained at least 30 million reads that aligned to 
RefSeq transcripts. Reads per kb per million reads (RPKMs) were 
calculated for each transcript by multiplying the number of aligned 
reads times 1,000,000 and dividing by the number of aligned reads 
for the replicate times the size of the transcript (in kb). 
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We then analyzed RPKMs to find significantly changing genes 
as previously described (Reddy et al. 2009). Briefly to account for 
the fact that lower-expressed genes are more variable between 
replicates, we used variance stabilization to normalize RPKMs. 
After variance stabilization, the normalized log ratios fit a normal 
distribution, and we assigned P-values accordingly. Once P-values 
were calculated for every transcript, the R package q value (Storey 
and Tibshirani 2003) was used to calculate false discovery rates, 
and we used a 5% cutoff to determine significantly changing 
genes. We performed this analysis comparing DMSO RPKMs to E2, 
GEN, or BPA in each cell line. To look for significantly enriched GO 
annotations, we used giprofiler (Reimand et al. 2007). 

Motif discovery and annotation 

To discover prevalent sequence motifs within the ESR1 bind- 
ing sites, we searched 50 bp on each side of the peak summit, 
which was reported by MACS (Zhang et al. 2008). We then used 
BioProspector (Liu et al. 2001) to perform de novo motif finding 
and look for sequence motifs with a width of 16 bp. The top five 
motifs for each ChlP-seq data set all matched the reported se- 
quence motif in Figure IB. Each peak was annotated with prox- 
imity to RefSeq transcripts. 

Quantitative PCR 

We performed qPCR to determine if long-term exposure to BPA or 
E2 influenced gene expression responses. To create cDNA, we used 
a modified version of the SYBR Green Cells-to-CT kit (Life Tech- 
nologies). We first washed each well twice with PBS, then added 90 
uX lysis buffer with DNase I, and incubated at room temperature 
for 5 min. Five jjlL stop solution was then added and incubated at 
room temperature for 2 min. To perform reverse transcription, 22.5 
juuL lysate was added to 25 jjlL RT mix and 2.5 uX RT enzyme and 
incubated at 3 7°C for 1 h and 95°C for 5 min. For qPCR, 3 uX of the 
RT reaction was added to 1 jjlL each of 10 |jlM primers (see Sup- 
plemental Table S5 for primer sequences), 10 uX PCR master mix, 
and 5 jjlL water. Biological duplicates were performed for each ex- 
periment, and expression measurements were normalized to CTCF 
expression levels. 

Data access 

ChlP-seq data can be accessed through the UCSC Genome 
Browser's site "Transcription Factor Binding Sites by ChlP-seq from 
ENCODE/HAIB" (http://genome.ucsc.edu/cgi-bin/hgTrackUi?db= 
hgl9&g=wgEncodeHaibTfbs) and through the NCBI Gene Expres- 
sion Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/) under 
accession number GSE32465. RNA-seq data can be accessed through 
the NCBI Gene Expression Omnibus (GEO) under accession number 
GSE38234. 
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